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Abstract 

We define a set of orthogonal functions on the complex projective space 
CP N_1 , and compute their Clebsch-Gordan coefficients as well as a large 
class of 6— j symbols. We also provide all the needed formulae for the gen- 
eration of high-temperature expansions for £7(iV)-invariant spin models defined 
on CP N - 1 . 
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1 Introduction 



Spin systems have been the subject of intense research for a long time. Indeed they 
describe the universal features of many different phenomena. For instance, the TV- 
vector model — a generalization of the Ising model — describes the critical behavior of 
dilute polymers, ferromagnets, binary fluids, liquid helium, and of the Higgs sector of 
the Standard Model at finite temperature [1, 2, 3]. From the point of view of quantum 
field theory, the iV-vector model provides the simplest example for the realization 
of a non-Abelian global symmetry. In particular, its two-dimensional version has 
been extensively studied because it shares with four-dimensional gauge theories the 
property of being asymptotically free in the weak-coupling perturbative expansion 
[4,5,6,7]. 

A generalization of the iV-vector model is provided by CP Ar_1 -models [8, 9, 10, 
11, 12, 13] (for recent work see [14, 15, 16] and references therein). Their two- 
dimensional version, besides being asymptotically free and invariant under global 
£7(./V")-transformations, shows some additional features that are expected to hold 
(though very difficult to prove) in the QCD case. Indeed, CP Ar " 1 -models have a 
local ?7(l)-gauge invariance, dynamical appearance of a linear confining potential 
between non-gauge-invariant states, and a non-trivial topological structure with in- 
stantons, anomalies and #-vacua. Besides being of interest in high-energy physics, the 
CP Ar_1 -models have applications in condensed-matter physics. Indeed, they describe 
models that have a complex vector order parameter (see e.g. [17]). 

Spin systems can be studied in a variety of ways: predictions can be obtained 
using field-theory methods — based essentially on perturbation theory [3] - or us- 
ing numerical techniques, e.g. Monte Carlo simulations or extrapolations of high- 
temperature series. The latter method (see for instance [18]) has proved very powerful 
in providing extremely accurate estimates of critical parameters. For spin systems, 
such as the iV-vector model, precise results [19, 20] are obtained even in the case of 
asymptotically free models, for which f3 c = +oo . 

Various methods have been used to generate high-temperature expansions for the 
N- vector model (for a review, see [21, 22]). The most straightforward method is the 
cumulant approach, or the closely related finite-cluster approach [23, 24, 25]. An 
improved method is the star-graph expansion, in which only star graphs need to be 
generated [26]. A further improvement was originally introduced by Joyce [27]. He 
proposes to expand the interaction term in hyperspherical harmonics. This allows 
a fast computation of the partition function for all star graphs. The main difficulty 
of the method is the computation of the Clebsch-Gordan coefficients, 6-j symbols 
and higher-order group invariants. Domb and coworkers [28, 29, 30] developed an 
algebraic method to compute the group factors associated to each graph. In Ref. 
[31] explicit expressions for the Clebsch-Gordan and the 6-j symbols were presented. 
Making use of simple algebraic rules [32] these general results allow the determination 
of all the group factors that are needed in the generation of high-temperature series of 
present-day length. These techniques have been recently used to generate long series 
in two and three dimensions [19, 32, 33]. We should mention that there exists another 
method that is extremely efficient in the generation of high-temperature series, the 
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linked-cluster technique [34, 35, 36, 37], which is an expansion in terms of free graphs. 

In this paper we will generalize the results of Ref. [31] to CP Ar_1 -models. We 
will use a representation of the hyperspherical harmonics 1 in terms of completely 
symmetric and traceless tensors. The computation of the Clebsch-Gordan coefficients 
is then a straightforward combinatoric exercise. Moreover, using the same technique, 
we are able to compute a very large class of 6-j symbols, essentially all symbols that 
are needed in the generation of high-temperature series of present-day length. 

The paper is organized as follows: in Sect. 2 we introduce the hyperspherical 
harmonics on CP N ~ 1 and discuss some general properties. General results and for- 
mulae for Clebsch-Gordan coefficients and 6-j symbols are reported in Sect. 3. In 
Sect. 4 we discuss the applications and report the high-temperature expansions of the 
Gibbs weight for the Hamiltonians that are used in studies of CP N ~ l models. The 
derivations of the results are reported in the Appendix. 

2 CP^^-Hyperspherical Harmonics 

In this section we introduce a class of orthogonal functions, which we call hyper- 
spherical harmonics, defined on CP N ~ l . These functions will provide the basis for 
expanding the Gibbs weight for U (iV)-invariant spin models. 

Let us first define CP^ -1 . Its usual definition is the following: consider w, z G 
C N — {0} and introduce the following equivalence relation: 

[Rell ] w ~ z if and only if there exists a G C such that w = az. 

Then 

C pN-l = ( C N _ { Q yy _ _ (2 J) 

A second definition is the following: consider vectors w belonging to the (2iV — 1)- 
dimensional complex sphere, i.e. w G C N such that W ■ w = 1, and introduce the 
equivalence relation: 

[Rel2 ] w ~ z if and only if there exists a G C, \a\ — 1, such that w = az. 
Then 

CP*- 1 = {z G C N : z ■ z = 1}/ ~ . (2.2) 

This second definition is obviously equivalent to the previous one. In our work we 
will use this representation. 

We mention a third definition that immediately shows that CP N ~ l is a symmetric 
space of rank one (see e.g. Refs. [39, 42]). Given p = (1,0, .. . ,0), notice that any 
z G C N such that z • z — 1 can be written as z — Up, with U G U(N). This 
representation is not unique: if U, V G U (N) are such that z = Up = Vp, then we 

1 These functions generalize the usual spherical harmonics for the three-dimensional rotation 
group [38]. Such a generalization is well-known for any symmetric space [39, 40, 41]. 
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can write p = U + V p . It is trivial to see that the set of R G U(N) such that p = Rp 
is isomorphic to U(N — 1). Therefore we obtain the result 

{z G C N : z-z = 1} ~ U(N)/U(N- 1) . (2.3) 

It follows easily 

CP N - 1 = SU(N)/S(U(N - 1) x 17(1)) . (2.4) 

It is interesting to observe that for N = 2, CP 1 ~ S 2 . The mapping is obtained as 
follows. Given z G C 2 such that ~z-z = \, consider 2 

= z a {aT p z^ (2.5) 

where a 1 are the Pauli matrices. Using the completeness relation 

{aT^fs = 2Stq-6ffi, (2.6) 

it is easy to verify that sV = 1, i.e. s G S 2 . Let us now show that the mapping is 
bijective. Given s G S 2 , if t = (0,0,1), there exists R G SO(3), such that R ij P = 
R i3 = s i_ Now given R e S0{3), it is well known that there exists U G SU{2) 
such that {U+^UYp = R ij (a j ) a p. Let us consider z a = U a p x p , with x = (1,0). 
Substituting this expression for z in Eq. (2.5), it is easy to see that one reobtains 
s* = R l3 . Finally, suppose that z x and z 2 are unimodular and both satisfy Eq. (2.5) 
with the same s l . Using Eq. (2.6) it is then easy to verify that 

|*i • z 2 \ 2 = i(l + sV) = 1 = |zi| 2 |z 2 | 2 . (2.7) 

This implies z 1 = e ia z 2 , for some a G 5R. Thus, Zi and z 2 are equivalent in CP 1 " 1 " 1 . 

Let us now define functions on CP N ~ 1 . It is obvious that a function / on CP N ~ l 
(defined in Eq. (2.2)) extends naturally to a function / on the complex sphere. Indeed, 
it is enough to define f(z) = f({z}), where {z} is the equivalence class (defined by 
[Rel2]) to which z belongs. On the other hand, a given function / on the complex 
sphere defines a function on CP N ^ if and only if f(z) = f(w) for z ~ w, or - 
in physicists' language — if f(z) is invariant under the [/(l)-gauge transformations 
z — > e ta z, i.e. 

f(e ia z,e- ia z) = f(z,z), (2.8) 

for any a. In the following we will identify the space of functions on CP N ~ l with the 
space of L r (l)-gauge-invariant functions defined on the complex (2iV — l)-dimensional 
sphere. We also introduce a normalized U (iV)-invariant measure 

dtt{z, z) = —— dzdz 5(z ■ z - 1) , (2.9) 

S 2 N 



where S 2 n is the surface area of the sphere S 2N 1 

(N-l)\ 



>2Af 



(2.10) 



In this paper we will use the usual summation convention over repeated indices. 
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Finally, we consider the unitary representation T(R) of SU(N) on L 2 (CP Ar_1 ) defined 
by (T(R)f)(z) = f(R~ 1 z). We want to find an orthogonal Hilbert-space decompo- 
sition of L 2 (C P N ~ l ) into subspaces such that the representation T(R) restricted to 
each subspace is irreducible. The needed decomposition 3 turns out to be precisely 
the decomposition of L 2 (CP N ~ r ) into eigenspaces of the Laplace- Beltrami operator 
C = C C p n - 1 - 4 I n f&ct, it can be proved 5 that: 

(a) The eigenvalues 6 of C are 

\ N ,k = k(N + k-l) > , (2.11) 
where k — 0,1,2, ... . The corresponding eigenspace En,k has dimension 7 

\N + k-2)\ l2 



M N , k = dimE N , k = (N-l)(N-l + 2k) 
and can be given several equivalent descriptions: 



k\(N- 1)! 



(2.12) 



(i) E Njk consists of the restrictions to the complex sphere z-z — 1 of the U(l)- 
gauge invariant polynomials P k (z,z) that are homogeneous of degree k in 
z,z E C N and that satisfy Laplace's equation ddP k = in C N . 

(ii) Consider the tensors Y Nk ^'"^{z) of rank k in z and z that are com- 
pletely symmetric under exchange of the a and j3 indices and that are 
traceless, i.e. vanish under any contraction of an a-index with a /5-index: 
YN,ktfpl'.'.'.p£ ( z ) = • We call these tensors spin-A; hyperspherical harmonics 
(they are described in more detail below). E^ tk is then spanned by the 
tensors Y Nk ^\"'^ as the indices range over the N 2k allowable values. 

(b) Each eigenspace E N ^ k is left invariant by T(R). Moreover, the representation 
T{R)\E N , k of SU (N) is irreducible. 

oo 

(c) L 2 (CP ) = E Nk (orthogonal Hilbert space decomposition). 

k=0 

To make all this concrete, we can write: 

YnaZa (*) = ^,k [z ai • • • z a *z Pl ...z h - Traces] , (2.13) 

where z E C N with z • z — 1, "Traces" is such that Y N>k is completely symmetric and 
traceless and 



1 

^N,k 



k\ 



(N-iy. ■ l ^ i4j 



3 See e.g. Refs. [39, 40, 41, 43]. 

4 The Laplace-Beltrami operator on CP 1 ^" 1 can be denned as the Casimir operator of the unitary 
representation of SU(N) on L 2 (CP N - 1 ). See e.g. Ref. [43]. 

5 See [41], Theorem 13.1, p. 231 and Corollary 14.4, p. 243. 

6 Note that our £ is the negative of the usual Laplacian, i.e. it is a positive-semide&nite operator. 
7 For a proof see Appendix A.l below. Note that A/jv,o = 1 an d A/jv,i = N 2 — 1 for all N. 
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Explicit examples are [the general formula is given in equation (A. 12)] : 

Y N , (z) = 1; (2.15) 

Y N ^(z) = v /iV(iV + 1)^-1^); (2.16) 



Y N , 2 ${z) = ly/N(N + l)(N + 2)(N + 3) 



2 

r a 8- — 



X 



: ' ZjZs — — — - (S^z^Zfi + 3 permutations) 



+ (j v + 1) 1 (jv + 2) (^ + ^) 



(2.17) 



We note that for iV = 2 we have X 2 ,k = + 1) and M 2 ,k — 2k + 1, as expected 
since CP 1 ~ S* 2 . Moreover, using the mapping (2.5), one can show that the Y's are 
linear combinations of the usual spherical harmonics. They can also be related to the 
real-valued hyperspherical harmonics (Inhere) N' k %k introduced in Ref. [31] by 

(Y sphcvc y 3 ) k ^(s) = 2-*/ 2 r 2j ^(; Z )K)\ i . . . {o*)\ h , (2.18) 

where s is related to z by Eq. (2.5) and we have used multiindex notation (a) k = 
aii . . . a k . For iV = 2 we can also prove the identity 

Y 2 *$:(z) = (-l) k Y 2 ^ k (w), (2.19) 

where w is defined by: w a = e a/3 Z/3, w a = e^z 13 . Here e a p = e a/3 is the completely 
antisymmetric tensor satisfying e 12 = e± 2 = 1. In group-theoretical terms, Formula 
(2.19) is related to the fact that S'[/(2)-representations are self-conjugate. 

The normalization is chosen so that the following orthogonality relation holds 
(see Appendix A.l): 

/ <m(z,z) Y N ^ k (z) Y N ^(z) = 8 kl I N , k $X , (2.20) 

where I^,k is the "identity" in the space of spin-/c symmetric and traceless tensors, i.e. 
the unique orthogonal projector onto the space of completely symmetric and traceless 
tensors of rank k. The trace of this operator is given by 

W(t ; S: = = dim E »* ' ( 2 - 21 ) 

as of course it must be. We remark that Y Nik (z) ■ Y N<k (z) = Y Nk ^ k k {z)Y Nk [a) k k ( z ) 
is independent of z [by U (N) invariance] , and hence 

Y N , k (z)-Y Ntk (z) = N N , k (2.22) 

by (2.20) and (2.21). 
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As stated in the theorem given at the beginning of this section, the hyperspherical 
harmonics are a complete set of functions on L 2 (CP N ~ r ). Thus any U(l)-gauge- 
invariant function f(z) can be expanded as 



/(*) = E Ut^j2:(*), (2-23) 

k=0 

where 



E Y N ^(z)Y Ntk ^ h (w) = 6(z,w) (2.25) 



For smooth functions this expansion converges very fast. Indeed, if /(z) is infinitely 
differentiable, then as k — > oo the coefficients of the expansion go to zero faster than 
any inverse power of k (see Appendix A.2). 

The completeness of the hyperspherical harmonics can be expressed through the 
relation 8 

oo 

1 N,k 

k=0 

where the 5-function is defined with respect to the measure dQ(z,z). 

Finally, let us consider a f/(Af )-invariant function / defined on CP 1 *' 1 . Let / be 
a function of two "spins" z, w, i.e. a function of \z • w\. We want now to compute its 
expansion in terms of hyperspherical harmonics. Using Schur's lemma we can write 

oo 

f(\z-w\) = ]T F N>k Y N , k (z) ■ Y N , k (w) . (2.26) 

fc=0 

The "Traces" terms can be dropped from either one of the Y's in the scalar product 
above, since the other Y is traceless. Also, since the scalar product is [/(AQ-invariant, 
we can rotate z to p — (1, 0, . . . , 0), and correspondingly rotate w to some v with 
\z ■ w\ — |p • v\ — Iv 1 ]. In this way we obtain 

Y N , k (z) ■ Y Ntk (w) = Y N , k (p) ■ Y N:k (v) = pL N , k Y N ^ i: i(v) . (2.27) 

Now Y N k \-\(v) can be expressed in terms of Jacobi polynomials 9 as (see Appendix 
AT) 



Y ^(v) - N ' k k { 1 1 - tl (2 2H) 



and therefore 



p r m e>\z-™\ 2 - 1) 

pf - 2 ,o )(1) 



Y N , k (z) ■ Y N , k (w) = N N , k k J N U,o )r : " • (2-29) 



8 Note that the normalization here follows directly from the one defined for (2.20). 

9 This expression corresponds to the relation between Yj and Legendre polynomials for the usual 
spherical harmonics. For definitions and properties for Jacobi polynomials see [44], pp. 1035-1037. 
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If we now invert Eq. (2.26) using the orthogonality relation (2.20) and Eq. (2.22) we 

F N , k = fdn(z,z) f(\z.w\) YN ' kiZ l YN ' k{w) . (2.30) 
This can be rewritten, using the [/(iV)-invariance of the measure and Eq. (2.29), as 

p(*-2,0)/ 2 . 1.2 _ i\ 

F N , k = j dQ(z,z) /(^|) k , N \ 2 \ 0) ] • (2.31) 

Now the integrand depends only on z 1 and we can integrate out the other coordinates. 
We finally obtain 

i p (N-2,o) (2 ^ 

F N , k — (N — 1) dt (1 — t)^ 2 /(Vf) fc p( „_ 2 V ~ J • (2.32) 

From the general properties of the hyperspherical harmonics we can derive the fol- 
lowing properties of the coefficients F N>k (for the proofs of properties 1 and 2, see 
Appendix A. 2): 

1. If /(*) is positive 10 for t G [0, 1], then \F Njk \ < F Nfi for all k ^ 0. 

2. If f(Vi) is smooth (i.e. C°°), then lim k n F Njk = for every n. 

fc^oo 

3. If f(t) = t 21 , then the integral in (2.32) can be performed explicitly and the 
coefficients F N>k are given by 11 

(N — 1)! (Z!) 2 



= { (N-l + l + k)\(l-k)\ ' " " ' (2.33) 

otherwise 
and are, in particular, always nonnegative. It immediately follows that for a 
generic function of the form 

oo 

/(*) = E fit*' (2-34) 

1=0 



if k < I 



the coefficients F^ tk are given by 



F N , k = £ /, F S • ( 2 - 35 ) 
Therefore, if all the coefficients j\ are nonnegative, then so are the F^ tk . 



10 More precisely, it suffices that / be nonnegative and not almost-everywhere-vanishing. 

11 See Eq. 11, p. 583, vol. 2 of Rcf. [45]. This result can also be obtained using Formula 7.391.3 
of [44], once one corrects the following misprint: T(a + 1) in the numerator should be replaced by 
I> + a + l). 
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3 Clebsch-Gordan and 6- j Symbols 

Let us now compute the Clebsch-Gordan coefficients for the decomposition into irre- 
ducible representations of the symmetric product of two generic representations. In 
general we can write 

V ^ k (-,\ v — ST r Wfc;(7);;(C)m v (e) m i \ /o i \ 

1 N,k,(fi) k \ Z ) 1 N,l,{S) t \ Z ) ~ L JV;fc,(, m ,(0) fc ;W i; ( £ ) m %,m,(() m l Z i ' ^.i; 

m 

Using the orthogonality relations (2.20) we obtain 

n (a)fc;(7)i;(f)m _ f tq/ -r\ v (a)k( \ v ( \ v (C)m / \ (n n\ 

U N;k,l,m,(f3)k;(S)i;(e) m ~ J a ^{ Z , Z ) 1 N,k,W)k ^ 1 N,l,(S) t \ Z ) 1 N,m,(e) m \ Z ) • \*-*) 

This integral can be computed explicitly. We get (see Appendix A. 4) 

(a)*;(7)li(C)m _ ( N ~ X ) ! VN,k VN,l VN,m k\ l\ ml 



r 

N;k,l,m,(/3) k ;(S)i;(e) m _ X _)_ ^ _|_ / _|_ m )1 

,, ,, ^T. 



'-ma x 



^ / k\ fl\ (m\ (a) k ;(p,)i(u) k _ z j (i)i;{r)j{&) k - h j (0™;(p)h(v)l-i (o o\ 



nun 



if \l — k\ < m < l+k ; the Clebsch-Gordan coefficient vanishes otherwise. In Eq. (3.3) 
h = m — l+i, j = m — k + i, i max = mm(k, I, l + k — m) and i min = max(0, l — m, k — m). 
In the following we will be interested in the scalar quantity 

C N; k,l,m = C N; k,l,m " ^N; k,l,m ■ (3-4) 

The general formula is reported in Appendix A. 4 [see (A. 60)]. A particularly simple 
case is m = I + k: 

r 2 - \r flN ' k fiN ' 1 a ^ 

N; k,l,l+k — J VN,l+k 2 ' y 6 - ) 

^N,l+k 

which can be obtained directly from (3.3), using (2.21). If k — 1 this gives 

c 2 N{N2 ~ l) ( N + l - 1 ) 2 (36) 

L N ; 1,1,1+1 ~ {N + 2 l) { I ) ■ W 

Notice that for N = 2, because of the symmetry (2.19), the Clebsch-Gordan coeffi- 
cients C N . k l m $)l'(])}'(e)™ — an d therefore also the scalar invariants C%. klm — vanish 
if k + / + m is odd. 

Let us now derive two important properties of the Clebsch-Gordan coefficients. 
Using their definition in terms of hyperspherical harmonics and the completeness 
relation (2.25) we can easily prove the crossing relation 

oo oo 

N;p,k,l,(S) p ;(e) k ;(C)i N;p,m,n,(a) p ;(a) m ;(T) n N;p,k,m,(8) p ;(e) k ;(a) m N;p,l,n,(a) p ;(Qi;(T)n ' 

p=0 p=0 

(3.7) 
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A second relation, which follows immediately from Schur's lemma, is 

n (a)fc;(/%(<0m n {l)k\{S)l\{v)n _ r J r (e)m;(/i)m /?2 /o o\ 

U iVifc,i,m,( 7 )fc;(«)i;(C)m ^^.l.n.laJklWcM, ~ °™» */■ J JV,m,(C) m ;H m °JV;fc,i,m ■ \*-°) 

Finally, using the completeness relation (2.25) and Eqs. (3.2), (2.22), it is easy to 
verify the identity 

oo 

E C%. ktltm = NniMn^. (3.9) 

k=0 

Let us now introduce the 6-j symbols (also called Racah symbols). They are 
U (iV)-scalars defined by 12 

1~> ( n L _. J _ f\ _ n (a)o;(/3)d!(7)c n Wal(5) 6 ;(e)e 

K, N {Cl,0,C,a,e, J ) - ^ N -a,d, c ,{rf) a ;(,\) d ;{lJ)c U Af;a,b,e,(a)a;H 6 ;(p) e 

P Wd;Wt;(C)/ /i H6;(m)c;(o-)/ mim 

X C 'JV;d,e,/,(^) d ;(6)e;( < 7) / C -JV;6, C ,/,(5) i ,;( 7 ) c ;(C) / ' ^- iU ^ 

The tetrahedral symmetry, which is enjoyed by the standard 6-j symbols [38], is 
trivially true also for our definition. 

We have not been able to compute a general formula for the 6-j symbols, but we 
have computed a very large class of special cases. This class is sufficient for computing 
the high-temperature expansion of the L r (A^)-invariant a-models, in general dimension 
d, up to rather high order. 

We begin by deriving a completeness relation for the 6-j symbols. Using Eqs. 
(3.7) and (3.8), and the definition (3.10) we get 

Y,K N (a,b,c;d,e,f) = -jj—C 2 Ndf C%. fif . (3.11) 

A general formula for IZjqil +p, I, p] r, m, k) for arbitrary /, p, r, m, k is reported in App. 
A. 4 [see Eq. (A. 70)]. It is particularly simple to derive all 6-j symbols in which one of 
the spins is 1. In this case we need to compute the coefficients TZn(1, l+a, k+b; k, I, m) 
with a,b — ±1,0. Using the tetrahedral symmetry one can immediately recognize 
that all cases can be computed using Eq. (A. 70) of App. A. 4, except when a = b = 0. 
In this case, however, we can use the completeness relation (3.11) to write 

1ZN(l,l,k;k,l,m) = —. C 2 N . klm C 2 N . lkk 

- /v iV,fc 

-K N (l,l-l,k;k,l,m) -K N (l,l+ l,k;k,l,m) . (3.12) 



4 Applications 

The formalism we have presented in the previous sections can be applied to the gen- 
eration of high-temperature series for U (iV)-invariant spin models defined on CP 1 ^^ 1 . 

12 In the terminology used in works dealing with high-temperature expansions [26] , the 6-j symbols 
are the group factors associated to the so-called a-topology. 
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The idea is to expand the Gibbs weight e~^ H for the given spin model in terms of 
hyperspherical harmonics, and then use their orthogonality to reduce the number of 
non- vanishing terms. We will consider here the Hamiltonians that are commonly used 
in the study of these models. 

The "quartic" Hamiltonian is given by 

H = -J2\z x -z y \ 2 , (4.1) 

(xy) 

where z x ■ z x = 1 and the sum is over all lattice links (xy). The Hamiltonian (4.1) is 
invariant under global U(N) transformations and under the local [/(l)-gauge trans- 
formations defined by 

z x - e ia * z x , (4.2) 

and therefore defines a theory on CP 1 ^" 1 . Notice that for N = 2, using Eq. (2.5), we 
have \z x ■ z y \ 2 — (1 + s x ■ s y )/2 . Thus the model with Hamiltonian (4.1) and N = 2 
is equivalent to the iV-vector model with N — 3. 

A second possibility is the so-called "gauge" Hamiltonian. In this case one asso- 
ciates to each lattice link a real gauge field 6 xy G [0, 2tc[ and considers 

H = -Re , £(z x -z y )e i0 ">, (4.3) 

(xy) 

where ~z x ■ z x = 1 and the sum is over all lattice links (xy). The Hamiltonian (4.3) is 
invariant under the local transformations 

z x ► e z x , 9 xy > 9 xy + ct x cn y , (4-4) 

and therefore defines a theory on CP 1 ^^ 1 . If one considers only correlations of the 
z-field, one can integrate out the gauge field 9 xy obtaining an effective Hamiltonian 

H = -W^g[I Q ^\z x .z y \)}. (4.5) 

' (xy) 

For the Hamiltonians (4.1) and (4.3), using (2.32)-(2.35), it is possible to compute the 
expansion coefficients F N;k . We need to compute F Njk for the functions exp (5 \z • w 
and I (P\z • w\). In the first case we obtain 

(N-l)\k\ Qk 



1 N,k 



(N-l + 2k)\ 



[3 k iF^k + I; N + 2k; (3) , (4.6) 



where 1F1 is the confluent (degenerate) hypergeometric function. 13 In the second case 
the integration gives 

F N , k = (N-l)\ ^ W-i(/3), (4.7) 
where I u is the modified Bessel function. 



13 



See [44], pp. 1058-1059. 
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A Properties of Hyperspherical Harmonics 



A.l Some basic formulae 

Let us begin by computing the dimension of the linear space E^^. This can be done 
by computing the dimension of the space of all completely symmetric tensors of rank 
k in two sets of indices, and then subtracting from it the number of independent 
trace conditions that have to be imposed in order to ensure the tracelessness of these 

tensors. The number of linearly independent symmetric tensors is given by ^ 
and the number of traces is given by (^^i 2 ) • Therefore we obtain 

Kf v t, [N + k-l\ 2 (N + k-2\ 2 .. t . 

N N , k = dimE N , k = u _ (A.l) 



(iV — l)(iV — l + 2Jfe) 



k - 1 

{N-2 + k)\ l2 
k\(N-l)\ 



(A.2) 



This proves Formula (2.12). For A^ = 2 we have N~2,k — 2k + 1, as expected on the 
basis of the isomorphism S 2 ~ CP 1 . 

Let us now compute the integral of monomials in z and z, i.e. of 

J cM(z,z)z a i...z ak zp 1 ...zp [ . (A.3) 

It is trivial to see that the integral is zero if k ^ I. To compute its value for k = I, 
we introduce an arbitrary complex vector A a , and define 

I k {A) = jdn(z,z)(A-z) k (z-A) k . (A.4) 

Since the measure dQ(z,z) is [/(AQ-invariant, we have Ik(A) = Ik{UA) for any 
U G U(N), so that Ik(A) depends only on A ■ A. Moreover, h(A) is manifestly a 
homogeneous function of degree k in A and A. Hence we must have Ik(A) = Jk(A-A) h 
for some constant Jk- Now, since z ■ z — 1, we have 

° ° h (A) = k 2 I k ^ (A) . (A.5) 



dA a dA a 

A recursion relation for Jk immediately follows: 



Using Jq = 1 we obtain the general solution 



/ - (A7) 
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Taking then k derivatives with respect to A and with respect to A in Eq. (A. 4) we 
obtain the well-known result 

/ d ^z, z)z^ . . . z^z Pl ...z Pk = ^ N _~f k)] -«+-) (A.8) 

where the terms in parentheses correspond to the k\ different pairings of the indices. 

Let us now prove the orthogonality relation (2.20). This is completely equivalent to 
proving that for arbitrary completely symmetric (in two sets of indices) and traceless 
tensors T^, k and Un,i we have 

J <M(z,z) (Y N , k ■ T N , k )(Y N ,t ■ U NJl ) = 8 kl T N , k ■ U N , k . (A.9) 

To see this, let us first use the definition (2.13) and let us notice that the "Traces" 
terms do not give any contribution due to the tracelessness of Tjv,jt and Um,i- Thus 
the l.h.s. in Eq. (A.9) becomes simply 

J dQ(z, z) z ai . . . z ak z^ . . . z^' z 7l . . . z lk z$ 1 . . . Zg t T N ^a) k ^jv,/,(/3) ( ■ 

(A.10) 



We can then use Eq. (A.8). The only non- vanishing contributions come from those 
terms which do not contain or Sg*; such terms exist only if I — k. In this last case 
there are (A;!) 2 equivalent contractions and we end up with 

(N-l)\(k\y st tt am 

Oki l^ N ,k ^ N z i + 2k)\ ' k ' ,h ~ ' k ' ,h ' y AAL ) 

Using the expression (2.14) for the normalization factor fiN,k, we obtain the orthogo- 
nality relation (2.20) for the Y's. 

The general formula for the hyperspherical harmonics [13] can be obtained by using 
the fact that they are completely symmetric and traceless. The complete symmetry, 
together with the needed transformation properties under U (N), implies an expansion 
of the form 

k 

Y N,k,{p) k ( Z ) = ^N,k A N,k;s P (k;s)(S)k ' ( A - 12 ) 

s=0 

where 

P {k ;smA z )= 8%...8%S* +1 ...*r i Zfi« 1 ...Zfi.+ permutations (A.13) 

and the number of permutations necessary to make P(k;s) completely symmetric is 
s\ (^j . Now we impose the tracelessness. We first note that 

= + (^^--iiwwt:^ > ( A - 14 ) 
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with the understanding that P( k - S ) = for s > k and s < 0. From the tracelessness 
of the hyperspherical harmonics we obtain the recursion relation 



A N , k;s + (N + 2k - s - 2) A N , k . a+1 = . (A. 15) 

Imposing the normalization A Ntk;0 = 1, we find 

(iV + 2A;- g -2)! 
^ = (iV + 2A:-2)! " (A ' 16) 

Let us now discuss the relation between the hyperspherical harmonics and the 
Jacobi polynomials. From Section 2 we know that Y N k \-\ (z) is the restriction to 
the complex unit sphere of a harmonic polynomial of degree k in z and ~z. Moreover, 
it depends only on l^l because of the gauge invariance z — > e ia z. Therefore it can 
be written as r 2k P k {\zi\ 2 /r 2 ) where r = \z\. Requiring the polynomial to satisfy 
Laplace's equation we get for P k (x) the equation 

d 2 P dP 

x(l - x) + (l - Nx) —± + k (N + k - 1) P k = 0. (A.17) 

dx 2 dx 

The regular solution of this equation is the Jacobi polynomial P% N 2,0 ' ) (2x — 1) (see 
[44, p. 1036]). The normalization is fixed by the requirement that 

^N,k,i.'.A ( z ) = ^N,k | z l | 2k + lower-order terms . (A. 18) 

We thus get (see Formula 8.962.1b in [44, p. 1036]) 

Y i- 1 ^ - ■ // kl (N-2 + k)\ (jv-2,o) / , |2 ,\ ( ,, q] 

which, using the fact that 

Pf-^d) = ( N ~ 2 k +k ) (A.20) 

gives (2.28). 

Note that we could have derived (A. 12) by using (2.28) and the expansion of the 
Jacobi polynomials. 



A. 2 Expansions in Terms of Hyperspherical Harmonics 

We want now to discuss the convergence of the expansion (2.23). We will begin by 

( a )k 



noticing that, given a generic complex tensor T N k ^ k , the hyperspherical harmonics 



satisfy the inequality 



T N , k ■ Y N , k (z) | 2 < (T Ntk ■ T* ;jfc ) N N , k , (A.21) 
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where the adjoint tensor T^ fc is given by 

M£"PW£. (A-22) 

and the overline indicates complex conjugation. Eq. (A. 21) follows immediately from 
Schwarz's inequality and (2.22). Moreover, equality in (A. 21) is possible only for 
those z's for which 

T N , k = 7 Y N , k (z) (A.23) 

for some constant 7. This requires first of all that T/v,fc be completely symmetric 
in both sets of indices and traceless. The constant 7 is easily obtained squaring the 
previous relation: 




(A.24) 



Now let us consider the special case T^, k = Y^ >k {x) with x = (1,0, . . . , 0) and 
k > 1. In this case we can rephrase Eq. (A. 21) in terms of Jacobi polynomials [using 
Eqs. (2.22), (2.28) and (2.27)] as 

Equality in this case is possible only if 

Y N , k (x) = lY N , k (z) , (A.26) 

with I7I = I. We will now prove that, for N > 2 this equation is satisfied only if 
7 = 1 and z = e l6 x. For N = 2 there is a second solution: 7 = (— l) k and z = e t9 x, 
with x = (0, 1). 

For b 7^ 1 we have from Eq. (A.26) 

V = Y N ,ktb(z) = VN, k (z 1 ) k (z b ) k , (A.27) 

that implies either z 1 = or z b = 0. Now, if z 1 = 0, from Eqs. (A.26) and (2.28), we 
obtain 

Pf- 2 '°>(1) = 7Pf- 2 ' 0) (-l). (A28) 

Now P fc (JV ~ 2 ' 0) (-r) = (-l) k , while ~ 2,0) (1) is given in Eq. (A.20). For N > 2 and 
I7I = 1, one immediately verifies that Eq. (A. 28) is never satisfied. Therefore z 1 7^ 0, 
so that z b = for all b ^ 1. The result immediately follows. For N = 2, using Eq. 
(2.19), it is easy to verify that both e l9 x and e l9 x satisfy Eq. (A.26). 

Thus, we have shown that equality in Eq. (A. 25) holds only for the case t — 1 
(respectively t = ±1) for N > 2 (respectively N = 2). We therefore have [since 
Pf- 2 ' 0) (1)>0] 

|Pf - 2 ' 0) (0l < Pf ~ 2 ' 0) (1) (A.29) 

for -1 < t < 1 and k > 1. 
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To discuss the convergence of the series (2.23) let us first notice that using the 
completeness relation (2.25) and Eq. (2.24) we get 



e k&hlS: = I dn ^*) i/(*)i 2 > ( A - 3 °) 

fe=0 J 

which is the Plancherel identity for harmonic analysis in CP^ -1 . Now let us consider, 
instead of /, the function C n f where C is the Laplace-Beltrami operator. In this case 
f k ^ is replaced by (\N,k) n fk > an d thus we obtain 

e (^) 2n k& 4(5: = / \£ n f( z )\ 2 ■ (A.3i) 

If now / is a C°° function, the integral is finite for all n. Thus the sum on the l.h.s. 
is converging for all n. Since Ajv,fc ~ k 2 for k — * oo, we get that k 2n fk • — > in the 
same limit, for every n. This implies that all coefficients fk decrease faster than any 
inverse power of k. To prove the convergence of the series (2.23) it is then enough to 
notice that \Y Njk (z)\ < (•A/jv,fc) 1//2 and that A/jv,jt behaves as k 2N ~ 3 for large k. 

Let us now discuss the properties of the coefficients F Njk in (2.32). The second 
property follows immediately from the previous discussion. We want now to prove 
that, if / (t) is positive for t G [0, 1], then \F N>k \ < F N;0 for k > 1. Indeed, from Eqs. 
(2.32) and (A.29) we get 



\F N , k \ < (N-1) fdt (l-tf- 2 f(Vt) 

J 



pf - 2 ,o )(1) 

< (N-1) f 1 dt (1 - tf- 2 f(Vt) = F Nfi ■ (A.32) 
Jo 

A. 3 Some Useful Contractions 

We are now going to report some useful formulae for the contractions of various 
products of hyperspherical harmonics. To simplify the results we will introduce some 
additional notation. 

Given three irreducible (completely symmetric and traceless) tensors T^,k, Un,i, 
and Vjv,m — respectively of rank k, I and m — we want to construct a scalar, i.e. a 
quantity that is invariant under U(N) transformations. It is easy to see that there 
are many ways of doing this, each one of them characterized by an integer %. We 
define 

[TN,k • U N j • Vjv, m ]j = 

— i+l — m(^)m+i — k T7 ( e )l — i (C)m — { A QQ\ 

1 N,k,(O m -i +l h)k- l+ i- m U N,l,(<*)i(e)i-i K A r ,™,(/3)fc- l (5) m+l -fc- l^'^J 

This contraction is defined for < i < (k, I) < m + % < k + /. In the previous 
expression, % counts the number of upper indices of T N;k that are contracted with the 
lower indices of U^^. Once i is given, all other contractions are completely fixed. 
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The definition (A. 33) can be generalized by writing 

rp (f)p TT T/ 

I N,k+p,(u) p ■ U N,l ■ V N,m\. = 

rp (H)p(a)i(f3)k-i TT {l)k-i + l-m{&)m + i-k T7 (Om-l+i ( A Q/H 

J Af,fc+p,H P (C) m -;+ ! (7)fc- l +i- m U JV,J,(a) i (e)j_ i l/ Af,m,(/3) fc _ l (5) m+ ,_ fe ■ l^-^J 

We also introduce a second useful notation: 

Z N*,$ h h (*)= z ai ---z ak z Pl ...z Pk . (A.35) 

Let us now compute various contractions that will be used in the determination 
of the Clebsch-Gordan and 6-j coefficients. 

(i) Bi(l, m; h, k). 

Let us define 

B 1 (l,m;h,k) = \Y n> i{z)-Y N ,m(z) ' ZN,l+m-k-h 

(z)] h , (A.36) 

for (l,m) > (h,k). By £/(iV)-invariance it is immediate to see that B\ does not 
depend on z. Using the property 14 

V I3{l)m-l( \ _ NN,mr L N,m-l v (j) m -l I \ (\ oq\ 

J v iV,m— 1 H-N,m 

we obtain for h > k 

B 1 (l,m;h,k) = X ":™ XN > 1 ^ h [Y N , h (z)-Y N , h (z)-Z NA _ k (z)} h . (A.40) 

Now we expand the first Yn,h{z) using Eq. (A. 12), and then use repeatedly Eq. (A. 39). 
We obtain finally (h > k) 

nn . Nn^ Mnjl i? Nth ^ (h\ fh-k\ 

^N, m Mm {N-2)\{N-2 + h + k)\ 

pLN,mVN,l (N-2 + h)\(N-2 + k)\- 1 • ' 

For h < k we can use the symmetry 

£i(Z,m;M) = #i(m,Z;M) . (A.42) 

Since the final expression, Eq. (A. 41), is symmetric in m and I and in h and k, it is 
valid also for h < k. 



14 This property is easily obtained recursively, by writing 



and noticing that 



Z N , m (z) ■ Y N , m (z) = = TT x k (A.38) 
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(ii) B 2 (s;i,l,m,n). 
Let us define 

B 2 (s;i,l,m,n)= [P(i ;s )(z) ■ Y NjTn (z) ■ Y N , n (z)\ . , (A.43) 

for < i < (to, I) < n + i < to + / and < s < I. 

The computation is simple. Indeed, using the definition (A. 13) of P(/ ;s )(z) we can 
easily express B 2 in terms of Si. If s < I — \m — n\ we obtain 

„ , . r ^ (i\ (n — m + i\ f I — A /to — n + I — i\ 

x r! (s — r)! I?i(to, n; n — I + i + s — r, r + m — i) , (A. 44) 

where r mir) , = max(0,i + s — /,i + s — l + n — to), r max = mm(s,i,i + n — to). Otherwise 
(if s > I — \m — n\) we get 

B 2 (s;i,l,m,n) = . (A.45) 

(iii) B 3 (i;l,m,n). 
Let us define 

£? 3 (*; Z, m, n) = \Y m (z) ■ Y N , m (z) ■ Y^z)^ , (A.46) 

for < i < (to, I) <n + i<m + l. 

Using Eq. (A. 12) we rewrite B 3 in terms of B 2 . It is easy to see that 

B 3 (i;l,m,n) = (A.47) 

if / < | to — ra | or I > to + n; otherwise we obtain 

l— \m— n\ 

B 3 (i;l,m,n) = fi N j A N j. s B 2 (s;i,l,m,n). (A.48) 

s=0 

Notice the following symmetries of B 3 : 

B 3 (i; I, to, n) = B 3 {1 — i; I, n, to) = B 3 (n + i — l;m,n,l) — B 3 (m + / — n — i; to, I , n) . 

(A.49) 

These symmetries are very useful to check the correctness of our result (A.48) since 
our method of computation gives an expression that does not have any of these 
symmetries in an obvious way. 

(iv) B 4 (i,j;p,r,m). 
Let us define 

B 4 (i,j;p,r,m) = J dfl(z,z) d£l(w : w) 

(z) ■ *V(«0]i (Y N , m (z) ■ Y N , m {w)) , (A.50) 
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for 0< <p< [(r + i),(r + j)}. 

Using the tracelessness of Y N ^ m {w) and Eq. (2.13) we can write 

B4(i,j;p,r,m) = /J, N , m [ dQ(z,z)dQ(w,w) 

(z) ■ Y N , r (w)] t (Z N , m (z) ■ Y N , m (w)) . (A.51) 

We can now perform the integral over z using Eq. (A. 8). It is easy to see that 
j',p, r, m) is non-vanishing only if |r — p\ < m < r — p + i + j. If this condition 
is satisfied we obtain (the integral over w can be dropped) 

fjL N>m (N -l)\i\ 3 \(m\y(r--p + 3 )\(r--p + i)\ 

BAi.y.p.r.m) = — ■ ; r- 

4V ' J '^' ; (N -1+r-p + i+j +m)\ 

_^ a\ b\ c! d\ e! f ! [Yn ' p{w) ' Y *' {w) ' Y ^)l +j -m- a - ( A -52) 

where b = j — a, c = r — p — m — a + i+ j, d = a + m — j , e=p + m + a — r — j, f = 
r+j-p-a, a min = max(0, j-m, r-p+j-m), a max = mm(j,r-p+j,r-p+i+j-m). 
We can now obtain the final result by using Eqs. (A. 46) and (A. 48). 

We note that, starting from the definition (A. 50) and using the completeness 
relation (2.25) and Eq. (A. 41), we obtain the sum rule 

r-p+i+j 

J2 B Sd;P,r,m) = B x {p,r\p-i,p- j) 

m=\r—p\ 

■N~N,r A/jv,p (N -2)\(N -2 + 2p-i-j)\ 



VN,rVN,p (N -2+p-i)\(N -2+p- j)\ 

(v) B 5 (i,l;p,k,r,m). 
Let us define 

B 5 (i,l;p,k,r,m) = j dVt(z,z) dVt(w,w)Y Nl f^ i (z) 

Ml 



(A.53) 



x 



y nmpM>) ' Yn ^ ' *W«0 L ■ Y N ,Jv>)) , (A.54) 



for < i < (p, k) < r + i < k + p. 

We expand Y Njk (z) using Eq. (A. 12) and notice that Y N j(z) can be replaced by 
(j,N t [Z N j(z). In this way B 5 is expressed in terms of £> 4 . We obtain finally 

B 5 (i,l;p,k,r,m) = n N ,k Vn,i [ ' 



, s — a , 



x (s — a)! a! Ajv,A:;s -B 4 (i — a + I, k — r + p — i — s + a + /; p + I, r, m) (A. 55) 

where a m j„ = max(0, s — k + i,s — k + i + r — p) and a max = min(i, s, r — p + i). 

Using the fact that B 4 (i, j;p, r, m) vanishes unless |r — p\ < m < r — p + i + j, we 
see that B 5 (i, l;p, k, r, m) — if m > k + I — s or m < \r — p — l\. 
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A. 4 Clebsch-Gordan Coefficients and 6- j Symbols 

Let us now discuss the computation of the Clebsch-Gordan coefficients (3.3). For 
arbitrary completely symmetric and traceless tensors Tjv,fc, Un,i and Vjv jm , we want 
to compute 

J dQ{z, z) ( T N , k ■ Y N>k (z) ) ( U N>1 ■ Y N>1 (z) ) ( V N , m ■ Y N , m (z) ) . (A.56) 
Using the definition (2.13) of the hyperspherical harmonics this reduces to 

m ' ^iV,m (■2-) ) • 

(A.57) 

We can now use Eq. (A. 8). For \l — m\ < k < I + m the integral becomes 

where h = m—l+i, j = m—k+i, i max = min(/c, /, l+k—m), i min = max(0, l—m, k—m) 
and we use the notation introduced in Eq. (A. 33). If k < \l — m\ or k > I + m the 
integral vanishes. Formula (3.3) immediately follows. 

We now discuss the computation of Cjj. fc z m for \l — m\ < k < I + m. Using (3.2) 
we rewrite 

— [ JO ( V -*\ r (a)fcl(7)i;(«)m V (ft/.\V (S)l f \ V (C)m/ \ 

^N;k,l,m - J OM\Z,Z) ^JV;fc,Z,m,(/3) fc ;(5)j;(C)m 1 N,k,{a) k \ Z ) 1 N,l,{l)A Z > 1 N,m,(e) m \ Z > " 

(A.59) 

Using (A. 58) we obtain 

^,2 (N — 1)\ /lN,k l~l>N,l ^N,m 

^N: k,Lm ~ 



(N — 1 + k + 1 + m)\ 



£ (£) Q (™)^!m! • YnM ■ > (A-60) 

where h, j,i min ,i max are defined as in Eq. (A. 58). The final result can be expressed 
in terms of B 3 (i] k, I, m) computed in Eq. (A. 48). 

We now describe another way of computing C%. klm . Using (3.2) we can write 

C%;k,l,m = J dtt(z,z)dtt(w,w) 

(Y N , k (z) ■ Y N , k (w)) (Y N , t (z) ■ Y N j(w)) (Y N , m (z) ■ Y N , m (w)) (AM) 

The integrand is only a function of \w • z\. Thus, using the U(N)- invariance of the 
measure, we can fix one of the two spins to an arbitrary value. Let us set w = x = 
(1,0, . . . , 0). We obtain, after integrating in dQ(w), 

C% k ,i, m = J dn(z,z) (Y N , k (z) ■ Y N , k (x)) {Y Ntl {z) ■ Y Ntl {x)) (Y N , m (z) ■ Y N , m (x)) 

(A.62) 
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and by using (2.29) we end up with 

C 2 N; k,l,m — (N — l)Af N , k X N ,lJ^N,m X 

y f 1 , f n .^-2 - 2 '° } (2t - 1) Pf- 2 ' 0) (2t - 1) PL N -^{2t - 1) 
A 1 } Pf- 2 ' 0) (1) Pf- 2 ' 0) (1) Pr- 2 ' 0) (1) ' 1 j 

Using the expansion of Pj; N ~ 2 ' \x) in powers of (1 — x) (cf. [44], p. 1036) we can write 
for \l — m\ < k < I + m 



C 2 N - k ,i, m = (N-l)Af N , k J^f N ,^ N 

(N - 

(N -2 + k)\(N -2 + s)\ 



* fk\(N-2)\(N-2 + k + s)\ l: l 



where 



,1 pW-2,0), x p (jV-2,0)( \ 

/(r, J, m) = 2— / <fc (1 - «)-*• % '* . (A.65) 

This integral can be evaluated explicitly (cf. Formula 2.22.18.2 of Ref. [45], vol. 2). 
For s < \l — m\ we have J(s; /, m) = 0, while for s > \l — m\ we obtain 



m 



(N -2 + l)\(N -2 + m)\ 

s + t)\ (iV-2 + t + m)!(iV-2 + t + s)! 



t =tl t!(m-t)!(s-Z + t)! (iV-2 + t)!(iV-l + t + s + Z)! 1 ; 

where t m i n = max(0, 1 — s). 

Finally, we report the computation of 7Zn{1 +P,P, r ', m , I, k), defined in Eq. (3.10), 
for \l + p — r\ < m < I + p + r, \l — k\ < m < I + k and \p — k\ < r < p + k. If 
one of these inequalities is not satisfied 1Zn(1 +p,p,r;m,l,k) =0, since CjV;s,t,u = if 
s > t + u or if s < \t — u\. 

Using Eq. (3.2), we rewrite 1Zn{1 +P,P, r; m, I, k) as 



J dQ(z, z) dVt{w, w) dVt(x, x) 



1 N,l+p,(f3) l+p \ W ) *N,l,(S)l \ Z ) r iV,p,(C) p W ^N;l+V,lMl+A~i)fMp 

x (Y N , m (w) ■ Y N , rn (z)) (Y N , k (x) ■ Y N , k (z)) (Y N , r (x) ■ Y N , r (w)) . (A.67) 
We can now use Eq. (3.3) to rewrite the above equation as 



/ dQ{z, z) dtt(w, w) dQ{x, x) \Y N)l+p (w) ■ Y N>l (z) ■ Yn^x)^ 
x (Y N , m (w) ■ Y N , m (z)) (Y N , k (x) ■ Y N , k (z)) (Y N , r (x) ■ Y N , r (w)) . (A.68) 
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Now we integrate over x using Eqs. (3.2) and (3.3). We obtain 

(N-l)\n 2 N n Nl n Nk n Nr r 

1 d\l(z, z) dll(w, w) 



(N- 1 +p + k + r)\fi NtP+ i 



X 



where h = p — r + i, j = p — k + i, i max = min(/c, r, fc + r — p) and i m j n = max(0, r — 
p,k-p). 

At this point we notice that the remaining term corresponds to B 5 , cf. Eq. (A. 54). 
We have therefore 

K„(l + vv r-m I k) - ( N ~ l )~ ^n, p Hn,i HN,k Hn,t 



(N — l+ p + k + r)\ Hn, p +i 



i- ma x 



x E ( Hnrjfc!Hp!B 5 (p-r + i;/;p,fc,r,m). (A.70) 



«=2 
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